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The widely used Double-Exchange model for manganites is shown to support various "striped" 
phases at filling fractions 1/n (n = 3, 4, 5...), in the previously unexplored regime of narrow 
bandwidth and small Jahn- Teller coupling. Working in two dimensions, our main result is that 
these stripes can be individually spin flipped without a physically relevant change in the energy, i.e., 
we find a large groundstate manifold with nearly degenerate energies. The two-dimensional spin 
system thus displays an unexpected dynamically generated dimensional reduction into decoupled 
one-dimensional stripes, even though the electronic states remain two-dimensional. Relations of our 
results with recent literature addressing compass models in quantum computing are discussed. 



PACS numbers: 75.25.Dk, 75.47.Lx, 75.25.-j, 71.10.Hf 

I. INTRODUCTION 

The manganite family of Mn-oxide materials has at- 
tracted considerable attention mainly due to the colossal 
magnetoresistance effect, where a magnetic field hugely 
enhances the electrical conduct ivity.l^ More recently, mul- 
tiferroicity has been found in some members of this fam- 
ily, particularly when chemical substitution reduces the 
bandwidth of the mobile electrons. Magnetic multifer- 
roicity, where a ferroelectric (FE) polarization is induced 
by the magnetic order, has lately been intensively in- 
vestigated.^^ Furthermore, manganites provide a fertile 
ground to study various ordered phases where magnetic, 
orbital, and charge degrees of freedom interact and com- 
pete. As a result of the large number of tendencies si- 
multaneously active on comparable energy scales, phases 
with very different ph ysic al properties can nevertheless 
be very close in eneigyW^ 

The competition of all these active degrees of freedom 
often leads to "striped" phases where electron motion is 
confined to one-dimensional (ID) subspaces. In mangan- 
ites at large hole doping, for example, electrons mainly 
occupy orbitals that point along either the x or the y 
direction. Spins then order ferromagnetically along this 
direction, because this favors the electronic kinetic energy 
via the double-exchange (DE) mechanism. In the other 
direction, where the electronic motion is suppressed, then 
DE is not active, and antiferromagnetic (AF) superex- 
change dominates. This establishes AF spin correlations, 
leading to the "spin striped" C-type AF phascfS^tini 
Other examples involve the ferromagnetic (FM) zigzag 
chains, which are antiferromagnetically coupled among 
themselves, that form the well-known CE-AF phase at 
half filling or the £'-AF phase in the undoped-limit.^'The 
ground state of these phases usually have collinear mag- 
netic order, i.e., alignment between any two spins is either 
perfectly FM or AF. Another aspect to note is that while 



the electronic kinetic energy is ID, the magnetic order is 
fully two dimensional (2D), with AF order between the 
chains mediated by superexchange (SE). In the absence 
of AF SE, a completely FM state with 2D electron mo- 
tion is energetically favorable. 

Recently, some of us predicted the existence of a 
new multiferroic phase, dubbed spin-orthogonal stripe 
(SOS) phase, located in the previously unexplored re- 
gion of quarter-hole-doping, small Jahn- Teller electron- 
lattice coupling, and narrow eg-electrons bandwidth.!^ 
In contrast to the "striped" phases mentioned above, 
the competition between FM DE and AF SE stabilizes 
non-coUinear magnetic order on some bonds, where the 
nearest-neighbor (NN) t2g spins are orthogonal to one 
another. As in an analogous half-doped noncollinear 
phase,'i2! this competition is expected to lead to a rel- 
atively high multiferroic critical temperature (~ 100 K). 
These phases can be visualized as composed of "thick" 
stripes, where all magnetic correlations within a stripe 
are collinear, i.e., AF or FM, while adjacent stripes dis- 
play non-collinear spins. 

In the present publication, we report that the SOS 
phase described above, as well as similar phases at 
dopings X = 1/n with integer n > 2, have another 
unexpected property, namely a very highly degener- 
ate ground-state manifold. We find that the spins of 
any collinear stripe can be rotated without a signifi- 
cant change in energy, as long as the spins in adja- 
cent stripes remain orthogonal. Spins in next-nearest- 
neighbor stripes can thus be at arbitrary angles relative 
to each other, implying that there is no magnetic order in 
the direction perpendicular to the stripes. Analogous di- 
mensional reduction effects on the magnetic order, where 
a higher-dimensional (2D in our case) spin system decou- 
ples into lower-dimensional (ID) subsystems, have been 
experimentally observed in other contexts such as near a 
quantum-critical point .'^Three-dimensional spin-ice sys- 
temJSl also show macroscopic degeneracy, and they can 
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decouple into 2D planes when a magn etic fi eld "switches 
oiT" some spins connecting the planes!^^'^ 

In contrast to spin ice, where a local symmetry al- 
lows two states for each tetrahedron, leading to a macro- 
scopic ground-state degeneracy and an extensive entropy 
proportional to the system size TV, the near degeneracy 
reported here is not quite macroscopic. The degener- 
acy involves flipping or rotating whole stripes, and the 
number of stripes in the 2D plane grows with the square 
root of the system size. This indicates that the emer- 
gent degeneracy of the ground state is intermediate be- 
tween local and global. Similar intermediate symmetry 
effects connected to dimensional reduction have been ex- 
tensively discussed in orbital-only models, as well as in 
the context of quantum computation. An example is the 
so-called "compass-model" where the x (y) components 
of the orbital pseudo-spins are coupled Ising-like in the 
X (y) direction. This model was originally introduced to 
capture the frustrated hoppings of the Cg orbitals of man- 
ganites.'^ Its low-e nergy states are ordered along only 
one dimensionj ^^ l ^^ l implying a symm etry t hat, again, is 
intermediate between local and global .I^^Eo] 

Even though the compass model was originally sug- 
gested to model the frustration inherent in degenerate Cg 
orbitals in 2D, the actual (partially frustrated) Hamil- 
tonian describing these orbitals does not exhibit such a 
high degeneracy. In fact, it has been shown that the pe- 
culiar properties of the compass model are sensitive to 
even rather slight modifications of the Hamiltonian and 
that a unique ground state with alternating orbital or- 
der sets in easily when the model is modified towards a 
more realistic description of 2D Cg systems.l^ In three 
dimensions, on the other hand, the model with a realis- 
tic Cg-orbital structure does decouple into planes^^^ as 
long as the magnetic order is fully FM. Once the mag- 
netic degree of freedom in included, however, the ground 
state turns out to show AF order along the z direction. 
This changes the relevant orbital Hamiltonian and as a 
consequence three-dimensional spin-orbital order with a 
non-degenerate ground state is stabilized. In addition to 
the eg orbital degrees of freedom in manganites, heavier 
elements with a strong relativistic spin-orbit interaction 
were discussed as a possible realization of compass mod- 
els.'^Sl However, here an isotropic Heisenberg term due to 
Hund's coupling would likely be present as well,'^ again 
inducing 2D order without 2^ degeneracy.'^ Finally, 
it was pointed out that the compass model describes 
properties desired for fault-tolerant qbit 8,1^21211 and cou- 
pled Josephson-junction arrays have been implemented 
for this purpose. 

The directional ordering described by the compass 
model is thus potentially very interesting, but dedicated 
Josephson-junction arrays appear so far the only systems 
showing such an effect. In the condensed-matter systems 
conjectured to display this physics based on particular 
features of their low-energy Hamiltonians, additional - 
even rather smalP^- perturbations, as often present in re- 
alistic descriptions of materials, have thus far been found 



to lift the high degeneracy. It has been suggested that 
the opposite route might work, i.e. that a more com- 
plex Hamiltonian, which does not itself have the appro- 
priate symmetry intermediate between local and global, 
might still s upport a ground-state manifold that has such 
symmetries .'22I32E1] mentioned above, some experi- 
mental evidence of such behavior exis t s,^^ and a similar 
effect is known in the case of spin-ice .^^^EH However, we 
are not aware of any model where "compass-like" behav- 
ior has been shown to emerge effectively for the ground- 
state manifold. We present here an unbiased numerical 
study showing this to happen in a Hamiltonian realistic 
for manganites, in a particular region of parameter space. 

Section[TT|contains the Model, which describes not only 
the eg orbital system at the origin of the compass-model, 
but also includes the spin and even the coupling to lat- 
tice distortions. In Sec. |III[ the numerical techniques are 
presented. In Sec. |IVt we show that the 2D spin system 
effectively decomposes into uncoupled ID stripes, while 
the electronic kinetic energy remains fully 2D. We also 
discuss orbital occupation and the relation between the 
magnetic degeneracy and dispersionless electronic states. 



II. MODEL 

The Hamiltonian considered here has been extensively 
studied in the past decade, and it has been shown to be 
very helpful to understand the complex behavior of man- 
ganites. In particular, the model reproduces the large va- 
riety of phases observed in manganites, e.g. A-type AF, 
C-type AF, FM, CE-AF, or £'-AF, and also the colossal- 
magnetoresistance regime.'^' The Hamiltonian is given by 
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The first term gives the kinetic energy of electrons in 
the two eg orbitals, containing the electronic hopping on 
NN bonds (ij) along the x and y directions (we study 
a 2D lattice for simphcity). The operator c}^ {c^^) cre- 
ates (annihilates) an electron at site i and at orbital a, 
where the orbital indices a and /? run over the dx2-y2 
and d^r^-z^ orbitals of the Mn ions. The orbital- and 
direction-dependent hopping parameters are given by 
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where the interorbital hoppings are negative (positive) 
on bonds along the the x (y) direction, and to = 0.2- 
0.5 eV defines the energy unit.'^l The Hund's coupling. 
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which links the itinerant electrons to localized t2g spins, 
is here taken to be infinite for simplicity, which implies 
that the Cg electrons' spin is always parallel to the local 
t2g spin. Neither Hund's rule coupling nor the electron 
spin thus appear explicitly in the Hamiltonian. This ap- 
proach leads to a modification of the bare hopping that 
is captured by a (complex) Berry phase factor 

n,, = cos I cos ^ + sin I sin ^ (3) 

^ 2 2 2 2 ' w 

which depends on the angles 9i, and 9j, 4'j that define 
the classical localized spins at sites i and Between 
sites with parallel (antiparallel) t2g spins, the Berry phase 
factor becomes one (zero), implying that the kinetic en- 
ergy favors parallel spins due to the DE mechanism. Be- 
tween non-collinear spins, i.e., spins with a relative angle 
between and 180 degrees, its absolute value is between 
and 1, and it can be negative or even complex. 

The second term of Eq. ([T]) describes the direct AF SE 
coupling between NN t2g spins. The third terms repre- 
sents the coupling of eg electrons with the lattice, via the 
Jahn- Teller (JT) {Q2 and Q3) and breathing (Qi) modes. 
A is a dimensionless lattice-electron coupling coefRcient. 
Only the x-y plane distortions are considered here, and 
if the overall lattice shape (square) can be assumed not 
to change, Qi can be set to —\/2Qz- The lattice normal 
modes Q2 and Q3 are related to shifts (5^^, ((5^ = 0) of 
the coordinates of the six oxygens surrounding each man- 
ganese via Q2 ^ ((5^ - <5'')/\/2 and Qg = -(P + 6^)1^/1. 
n-ia = c\^c^^ and rn = ni^^2_y2 + ni_3^2_r2 are density 
operators; while 

< = {ni,x^-y2 - nj_322_^2)/2, (4) 
'''i = ("^l^^-ySC- 3^2_r2 + cJ,3z2-r2'^i,x2-y2)/2, (5) 

denote the orbital pseudospin operators, similar to the 
Pauli matrices for spins.'^' The lattice distortions are thus 
coupled to the orbital degree of freedom. Undoped man- 
ganites correspond to a filling of one electron per Mn, 
while doping x and filling n are related via n = 1 — x. 

III. METHODS 

Markov-chain Monte Carlo (MCMC), zero-T opti- 
mization, and variational methods were used to study 
the ground-state and low-temperature properties of the 
Hamiltonian Eq. Q. This Hamiltonian couples non- 
interacting electrons to the classical variables ql = 
{6^ ,6f,S^} describing the localized spins and lattice 
distortions. For any set of (f s, Eq. ([!]) is diagonalized us- 
ing standard library routines, and the free energy of the 
electronic system can be easily evaluated via the usual 
equations from statistical physics. The energy is then 
used in a conventional MCMC procedure to determine 
whether a new set of -0*^ should be accepted, as detailed in 
various publications.'^' Since care must be taken that the 



results are independent from the initial configuration and 
that thermal equilibrium has been reached, such MCMC 
simulations are quite CPU-time consuming. Their huge 
advantage is that they are unbiased, meaning that, for 
long enough run time, they will converge to the true rel- 
evant state of the cluster under study, regardless of the 
initial state used. We employed this method on 6 x 6 
lattices to obtain the phase diagram in the A- Jaf plane, 
and on 12 x 12 lattices for a selected number of points. 

The MCMC is complemented by the zero-T optimiza- 
tion method where the (fs are optimized to reach the low- 
est possible energy by employing the derivatives 
as detailed in Ref. 1331 This optimization method is par- 
ticularly useful around a local energy minimum, where 
it reaches higher precision than the MCMC and is, thus, 
efficiently combined with MCMC, which is better at find- 
ing the global minimum. Finally, we also complement 
these studies by a variational comparison of the energies 
of fixed configurations of classical spins and lattice dis- 
tortions. While the variational approach is not unbiased, 
because only chosen configurations were combined, it is 
valuable, because far larger lattices can be reached; it 
was employed, e.g., to verify that various SOS1/3 config- 
urations indeed are practically degenerate. In addition 
to L X L squares, this approach was also employed on 
\/2Li X ^/2L2 rectangles. Periodic boundary conditions 
were used in MCMC and optimization, and for varia- 
tional energy comparison. To calculate selected observ- 
ables, such as the spin-structure factor and the density of 
states, we additionally used twisted boundary conditions 
leading to a denser k mesh, as described in, e.g., Ref. [Ml 

IV. RESULTS 

Motivated by the recent prediction of the x = 1/4 
SOS phase,lID and of a similar phase at half-dopinji^l -^^^rg 
have investigated other fillings x = 1/3, 1/5, . . . , focus- 
ing here on parameters relevant to narrow-band mangan- 
ites. The discussion below mainly focuses on results for 
X = 1/3, where the SOS phase has narrower stripes than 
for smaller x and can thus be investigated on smaller clus- 
ters. A variational energy comparison and some MCMC 
studies were also performed for n > 5, leading to analo- 
gous results. 

A. Phase diagram and highly degenerate ground 
state manifold 

Figure [ija) shows a snapshot obtained in a MCMC 
run for Jaf = 0.1 9to and A = 0, where the spins hap- 
pened to lie almost within the x-y plane. It illustrates 
the SOS-phase expected for a fiUing of 1/3, actually just 
one of its realizations, see the discussion below. As dis- 
cussed for X = 1/4,^^^^ the SOS phase consists of domains 
of the £'-AF phase of undoped manganites [illustrated in 
Fig.[l](b)], with spins rotated by 90° between neighboring 
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FIG. 1: (color online) (a) Monte-Carlo snapshot for Jaf = 
0.19to, A = 0, and /3to = 1200, in the regime of the SOS1/3 
phase. It can be visualized as composed of domains of the 
-E-AF phase illustrated in the cartoon (b). NN domains 
("stripes") are separated by dashed lines in (a), they have 
spins at right angles. In (b), dashed lines indicate the zig-zag 
FM chains of the iJ-AF phase; in (a) shading (color) indicates 
a few of the short segments of the i5-phase zigzag-chains that 
survive in the SOS1/3 phase, later called "arrows". 




FIG. 2: (color online) Zero temperature (T = 0) energies (per 
site) of several magnetic phases, at doping x = 1/3 and A = 0, 
vs. Jaf. Dots represent MCMC results at T = 0.002 on a 
6x6 lattice that closely follow the variational results. S6 
denotes a spiral phase with a period of 6 lattice spacings,'^ 
while a 12 X 12 lattices shows one with period 12; this phase 
thus converges to FM order with increasing lattice sizes.'^ 



domains. Each domain can be visualized as one "stripe" , 
see Fig. [ija) , and the spins are then coUinear within each 
of the stripes (regions between a pair of dashed lines), 
but between stripes they are orthogonal to each other. 
In Fig.[2j the groundstate energy of the SOS1/3 phase is 
compared to that of various other phases in the absence 
of electron-phonon coupling (A = 0), and it is clear that 
the SOSi/3 phase has the lowest energy for a range of 
Jaf- This is corroborated by data points obtained with 
unbiased MCMC, which closely follow the variational en- 
ergies, indicating that the true ground state has been 
found. 

The phase diagram including electron-phonon coupling 
A in addition to Jaf is given in Fig. |3j and it shows 



FIG. 3: (color online) Zero-T phase diagram of model Eq.(l) 
at a;=l/3, varying A and Jaf- Variational technique results 
are shown with shading (colors). The FM DE metallic phase 
dominates; large Jaf stabilizes AF phases. The "C1/3E2/3" 
phase is a varianlP^ of the CE phase at half filling.'^ The 
SOSi/3 phase at intermediate Jaf and small A is analogous 
to the a;=l/4 SOS1/4 phase and is highly degenerate. Full 
dots show where MCMC (for 6x6 sites) confirmed the results; 
open dots are MCMC results that remained inconclusive due 
to metastabilities caused by phase competition. For details 
see Ref . HH 



that the SOS 1/3 phase remains a stable ground state for 
A < 0.7. As it may be expected, DE drives a FM metallic 
state at smaller Jaf, while large Jaf stabilizes a fully G- 
type AF phase. At large A and intermediate Jaf, the 
SOSi/3 phase is replaced by a variant of the exotic but 
well studiecPSl CE-phase. 

However, our main result is that our calculations have 
revealed that this ground state is not unique: if all spins 
within one stripe are flipped by 180° (i.e. inverted), the 
total energy remains nearly unchanged. This is illus- 
trated in Fig. Qa-b) , where two almost degenerate vari- 
ants of the SOSi/3 phase are illustrated. Taking into ac- 
count non-coplanar spin configurations as well, the spins 
in each stripe can in fact be rotated by any angle, as long 
as the spins in adjacent stripes remain orthogonal. 

On square clusters, the difference in energy per site 
before and after a stripe spin-flip is merely ~ 10~^to for 
flipping one stripe of a 12 x 12 lattice. This is three orders 
of magnitude smaller than the energy differences with 
other states shown in Fig. [2j similar conclusions were 
reached using tilted rectangles. Finite electron-phonon 
coupling A > as well as finite Hund's rule coupling 
>^Hund < 00 leave the near degeneracy intact. Even rather 
small Hund's coupling 6io ~ 1-2 eV increases the energy 
split only to ~ 5 • 10~^io, still more than two orders of 
magnitude smaller than the energy differences to other 
phases (A > O.Ttg eventually drives a transition to the 
C1/3E2/3 state, see Fig.jsj) 
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FIG. 4: (color online) Cartoon for the spin patterns for two 
quasi-degenerate states of the SOS1/3 phase (a-b). Lines be- 
tween domains with orthogonal spins illustrate the periodic 
pattern of Berry phases: Thin lines denote positive sign, thick 
lines negative. The pattern in (b) is reached by flipping all 
spins in the top black stripe of (a), which entails changes in 
the signs of the Berry phases, but induces an energy differ- 
ence of only 10~^to. The dotted lines connecting three ferro- 
magnetically aligned spins indicate 'arrows' of the i5-phase, 
see Fig.[l] In the pattern in (c), 'arrows' in adjacent domains 
point in opposite directions. The energy per site in this phase 
is 10~^to higher than that of the SOS1/3 phase; flipping spins 



of one stripe induces energy differences of « 5 ■ 10 
least 50 times larger than between (a) and (b). 



to, i.e., at 
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FIG. 5: (color online) Spin-structure factor ^(k) along the 
line (tt, 0) to (0, tt), calculated for the SOS1/3 phase by aver- 
aging over all 2^® degenerate realizations with stripes along 
one direction of a 24\/2 x \/2 cluster (A=0). S(k) ^ for all 
k except the line running from (vr, 0) to (0,7r). Results at x 
— 1/4, 1/5, and 1/6 are also shown. 



B. Magnetic, orbital and charge patterns 

Effectively, the 2D system decouples into ID stripes, 
whose direction can be rotated independently from the 
other stripes as long as spins in NN stripes are at right an- 
gles. The relative orientation of stripes at larger distance 
is thus arbitrary and the spin structure factor S'(k) is fi- 
nite along a whole line in momentum space, see Fig.[5][b), 
similar to results for the compass model.'^ The 5(k) 
modulation is due to the width of a double stripe. 



(b) 



FIG. 6: (color online) (a) Zero-T spin and orbital ordering 
of the SOSi/3 phase, the ground state for 0.17 < Jap 0.21 
(A = 0), see Fig. p] For each site, the linear combination 
\a) = cosq:|3z^ — r) + shia\x^ — y^) is depicted that has the 
highest density, i.e., for the value of a maximizing rii^a- In 
(b), the orbital with maximal density is shown for the mod- 
ified SOSi/3 phase with the spin pattern given in Fig. [4|c). 
While the depicted orbital is the one with the highest density, 
these cartoons do not fully describe the orbital states, as some 
density is also found in the orthogonal orbitals, see text. 



Since spins on all bonds between stripes are orthog- 
onal, the absolute value of the hopping connecting the 
stripes is the same in all realizations of the SOS1/3 phase, 
with = l/\/2 obtained from (l3|. However, having 



the same \VLij\ is not enough to establish such a degener- 
acy, as the complex phase of the Berry phase in general 
cannot be neglected. In the "flux" phase,'^NN spins are 
always orthogonal and flipping a spin would not change 
this; but there is nevertheless a unique ground state sta- 
bilized by the Berry phase. Similarly, flipping a stripe 
of the "SOS1/2" phase at half dopingl-^ preserves the ab- 
solute value of all hoppings, yet costs far more energy 
than flipping a stripe of the SOS1/3 phase. In order to 
show that the groundstate degeneracy is not due solely to 
having orthogonal spins, it is illustrative to analyze the 
"modified" SOS1/3 phase depicted in Fig. |4|c). Like the 
actual SOSi/3 phase, this phase is made up of domains of 
the i?-AF phase that are orthogonal to each other. The 
only difference is that the "arrows" formed by the FM 
spins point in opposite directions in adjacent stripes in 
this modified phase. Even though the spins along the 
boundaries between stripes are the same in both phases, 
the modified phase has a higher energy per site by 10~^ioj 
which is larger that the energy difference to the C1/3E2/3 
phase, see Fig. [2j indicating that the internal composi- 
tion of the stripes matters as much as their boundaries. 
Moreover, the modified phase does not show such a near 
perfect degeneracy, as flipping a stripe on a 12 x 12 lattice 
costs ~ 5- IQ^Hq, and while this is not a large energy dif- 
ference, it is at least 50 times as much as for the SOS1/3 
phase. 

The dominant orbital occupancy and its relation to the 
spin pattern is shown in Fig. |6(a) The building block of 
the SOSi/3 phase is an arrow made of three FM spins, 
with a 'center' and two 'wings' that point in the x and y 
directions. For each site, one can calculate the dominant 



FIG. 7: (color online) (a) Density of states N{uj) (A=0) 
for two patterns of the SOS1/3 phase, the phase shown in 
Fig. ^c), and the SOS1/2 phase.'^ The arrow indicates the 
dispersionless states of the SOS1/3 phases. 



linear combination of the two Cg orbitals, i.e., the orbital 
with the highest density. At the wings of the arrows, 
these are the directional orbitals pointing to the center, 
because they can maxim ize the kinetic energy along the 
FM bond,'i'see Fig. 6(a) they are half filled. On the cen- 
tral site, the x'^ — orbital dominates due to its large 



(b) 



overlap with adjacent sites. However, it turns out that 
the electronic configuration of the SOS1/3 phase cannot 
be fully characterized in terms of singly occupied orbitals, 
which is reminiscent of the edge sites in the CE-phase.'^® 
Even though the depicted orbitals have the highest den- 
sity, some electronic weight is also found in the orthog- 
onal states and Fig. 6(a) only partially describes the or- 
bital state. In the center, the explicit occupancies are 
thus n^2_y2 — 0.4 and n2^2_^2 ~ 0.28, and in the 'wings' 
along xjy^ we find ny2_^2 j ^2_^2 = 0.24 in addition to 
n32.2_r2/3y2_^2 = 0.5. The total density is the sum of 
the two orbital densities on each site, and shows only 
weak stripe modulation, with n — 0.68 in the center and 
n = 0.66 in the wings. 

The relatively large density ny2_^2 = 0.24 in a 'wing' 
pointing along x may seem surprising, as the — or- 
bital cannot hop along the FM bond directed along x to 
the center.ii' However, it can hop with an only slightly 
reduced amplitude = l/-\/2 to an adjacent stripe. 

This process becomes even more important in the modi- 
fied SOSi/3 phase of Fig. |4|c), whose dominant orbitals 
are shown in Fig. 6(b) In this case, the 1 / \/2 inter-stripe 



hopping connects two directional 3a; 



(or 3y 



orbitals with a large overlap. We will discuss in the next 
Section |IV C| why the hoppings between stripes cannot 
be neglected in either phase. 



C. Two-dimensional kinetic energy and 
dispersionless states 

As discussed above, the effective hoppings with their 
Berry phases are not the same in the various (almost) de- 
generate states of the SOS1/3 phase. It may be tempting 
to believe that one can find a local gauge transforma- 
tion that transforms the sign patters of one realization 



2 
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1 
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FIG. 8: (color online) Spectral density A(k,(jj) for two pat- 
terns of the SOSi/3 phase at A = 0. The Lj; x Ly = 24 x 24 
lattice is used, and the E-KF domains run in the (1, 1) direc- 
tion, from (0,0) to (24,24). Periodic boundary conditions are 
employed, and peaks were broadened with a Lorentzian with 
a width 5 = 0.02fo. 



into that of another, with a flipped stripe. In this case, 
the electronic Hamiltonians would be equivalent and thus 
have the same eigenenergies, explaining the degeneracy. 
However, the density of states shown in Fig. [7] is very dif- 
ferent for different SOS1/3 patterns, indicating that the 
eigenenergies of their electronic Hamiltonians are in fact 
quite different. 

That the different SOS1/3 conflgurations have different 
electronic Hamiltonians can also be inferred from the one- 
particle spectral density v4(k, lo) shown in Fig. [sj One 
aspect to note is that the occupied states do not show 
any (quasi-) ID character, and the states are dispersive 
both along the stripes [(0, 0) - (tt, tt)] and perpendicular 
to them [(0,0) - (tt, — tt)]. In fact, the kinetic energy 
between the stripes is necessary for the SOS 1/3 phase, as 
the orthogonal spin arrangements are stabilized by the 
competition between the FM DE and AF SE processes. 
Electrons can not hop directly along the stripes, as NN 
spins are AF ordered, but they move along this direction 
via neighboring stripes, leading to a somewhat weaker 
dispersion than perpendicular to the stripes. 

The spectral density of the SOS1/3 phase in Fig. |8|can 
be compared to that of a toy model where all hoppings 
are replaced by their absolute values. As can be seen in 
Fig. [9] the spectral density then consists of six coherent 
bands (the dispersionless band at a; = is doubly de- 
generate) corresponding to six states formed by the two 
orbitals of three sites within each arrow. The results for 
the full model in Fig. [8] shows some remnants of these 
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(0,0) (0,Jt) (71,0) (71,71) (0,0) (7t,-7[) 



FIG. 9: (Color online) Spectral density A(k, uu) for an artifi- 
cial reference model where all hoppings are replaced by their 
absolute values; this phase has the same energy as the SOS1/3 
phase. Parameters are as in Fig. |8] 



bands, especially around (0,0) and (7r,7r). At other mo- 
menta, however, the bands split or the weight even ap- 
pears incoherent. This reflects the additional modulation 
of the hoppings by the complex Berry phase, which can 
lead to superstructures and to effective disorder. As one 
can infer from the density of states, see Fig. [7| spectral 
weight is not only shifted to different momenta in differ- 
ent SOSi/3 configurations, but also transferred between 
different energies. 

The dimensional reduction unveiled here appears 
linked^to the dispersionless edge-states of the E-AF 
phase,'^'2II which are unoccupied in the SOSi/^ phases 
[see the delta peak in the density of states, Fig.lTj and the 
dispersionless states in A(k, w) in Fig.js]. While the onc- 
particle energies of the occupied bands change, the dis- 
persionless states are unaffected by flipping a stripe. The 
"SOS1/2" phase,^^! and the modified phase of Fig.[4|c) do 
not have such dispersionless states and do also not have 
degenerate states, and we will thus analyze the disper- 
sionless states and their impact on the SOS1/3 phase. 

In finite E-AF clusters, the dispersionless states are 
z'^ — and 
edge. These 



z orbitals localized on sites along the 



orbitals are localized be- 



cause their orbital symmetry only permits them to hop 
to a single site, their neighbor in the a;-(t/-)direction, but 
this bond is AF, and hopping is thus suppressed by the 
magnetic order. In the bulk, the bond in the opposite 
direction would be FM, and the orbital can thus delo- 
calize, but this bond is missing at the cluster edge. If 
one considers a single stripe of the SOS1/3 phase as an 
isolated _E-AF domain, then the electrons can only delo- 
calize within (i) both orbitals of the central site and (ii) 
the directional orbitals pointing towards it on the two 
wings [these orbitals are shown in Fig. |6(a)| . The or- 
thogonal (planar) orbitals on the wings do not hybridize 
with any other orbital of the stripe and are the localized 
edge states, which are empty at x = 1/3. 
As discussed above in Sec. 



IV B however, hopping be- 



tween the stripes is not suppressed and the arrows on ad- 
jacent stripes are thus coupled. As illustrated in Fig. [TOj 
electrons in the planar orbitals can hop to the directional 




FIG. 10: (Color online) Two orbital building blocks of the 
SOSi/3 phase. Black arrows give the SOS1/3 spin pattern 
and dashed lines indicate the stripes. Orbitals shown with 
shading/lines illustrate the unit cells corresponding to two 
arrows in adjacent stripes. Site A is a 'center' site with occu- 
pancy in both the — and the 3^^ — r'^ orbitals, site B is a 
'wing' site, only the more highly occupied — orbital is 
shown. Site C is also a 'wing' site, but here both orbitals are 



shown, the planar y 



orbital is drawn with lines and the 



directional — with shading. An electron in the planar 
y^ — orbital can only hop along -\-y to site B of the adjacent 
stripe: It cannot hop along the a;-direction due to its orbital 
symmetry,'^' and it cannot hop along the —y bond due to the 
AF spins. 



orbital on the wing site of an adjacent stripe. This leads 
to the subsystem consisting of the six orbitals drawn with 
black lines in Fig.[TOj After taking into account reflection 
symmetry with respect to the central site, the even and 
odd subspace each yield a 3 x 3 matrix 



/o 
















_ V3to 
2\/2 





V^*0 a — i<t> 

2y^^ 






(6) 



where subscripts e/o denote the even and odd subspace, 
ie,o is given by ^ for the even and for the odd case, 
and (j) denotes the complex Berry phase. It can be easily 
seen that such a matrix always has one eigenvalue e = 
and that the corresponding eigenvector only lives on the 
first and third sites. Here, this implies that the direc- 
tional orbital on the wing sites has no overlap with dis- 
persionless states. The only process connecting the build- 
ing blocks described by Eq. (6| is the hopping — 
between two directional orbitals sitting on wing sites be- 
longing to adjacent stripes, e.g., between the 3?/^ — 
orbital on site B and the 'ix^ — orbital on site C in 
Fig. 10 Since it only acts on the directional orbitals. 
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which are not involved in the e = states, the latter 
are not afTected and remain dispersionless in the SOS1/3 
phase. In the "SOS1/2" phasef^^l the E'-phase domains 
are too narrow to support dispersionless states and in 
the modified phase of Fig. ^c), where the "arrows" of 
the E-phase point in opposite directions in adjacent do- 
mains, the different orbital pattern shown in Fig. |6(b)| 
also prevents similar dispersionless states. 

Such a protection of a large ground-state degeneracy 
by dispersionless states is reminiscent of spin ice, where 
flat bands enforce the ice rules."*" The notion that the 
protection of a degeneracy is caused by the dispersionless 
states at w = can be motivated by observing that the 
lack of dispersion implies that the states are completely 
localized and do not communicate with other, dispersive, 
states. The energy of the occupied subband centered 
around cj = — e « — 1.2io can then only be affected by 
processes involving the symmetric unoccupied states at 
oj = +e. Treating the hybridization between these two 
subbands in perturbation theory, valid for large energy 
differences 2e > the second order contribution can 

2v 2 

be shown to be independent of the hopping's complex 
phase and the third order contribution drops out entirely. 
If hybridization sensitive to the Berry phase then occurs 
almost exclusively within the occupied subband, the in- 
dividual one-particle energies seen in the density of states 
can be changed, but the total energy cannot, because all 
energy gained by one state is lost by another. 

V. DISCUSSION AND CONCLUSIONS 

The E-AF phase, which provides the stripe build- 
ing bloc ks, d isplays FE polarization due to exchange- 
striction.Sillll Additionally, the many bonds with non- 
coUinear spins leads to a sizable Dzyaloshinskii-Moriya 
interaction. We found that the FE polarization caused 
by the latter is larger in some configurations and smaller 
in others. This causes an additional energy difference. 



but as long as it is small compared to the energy sep- 
arating the SOSi/3 ground-state manifold from excited 
states, this should not qualitatively alter the physical be- 
havior that we described here. On the other hand, the 
different FE properties might provide a handle to ma- 
nipulate the stripes. It should finally be noted that the 
multiferroic properties of a large collection of nearly de- 
generate states differing by the FE polarization has not 
been investigated thus far and may be highly non-trivial 
due to potential interference effects. 

To conclude, we have studied a model Hamiltonian 
appropriate for narrow-band manganites at small Jahn- 
Teller coupling using unbiased numerical techniques. We 
found that the spins of the 2D system spontaneously un- 
dergo dimensional reduction into ID stripes for dopings 
1/3, 1/4, ... Adjacent stripes always have spins at a 90° 
angle, but the spins of stripes at larger distances are not 
correlated. The electronic kinetic energy, on the other 
hand, remains fully 2D. This indicates that the mangan- 
ite results described here induce the spins to behave in 
a way analogous to the (orbital pseudo-) spins in the 
so-called compass model, which is also used to describe 
protected qbits. However, in contrast to the compass 
model, our model Hamiltonian does not commute with 
the corresponding symmetry operators, and the effect is 
thus an emergent property of the groundstate manifold. 
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